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Abstract. We present a statistical approach for estimating the age of Globular Clus- 
ters by measurement of the likelihood between the observed cluster sequences in the 
Colour-Magnitude Diagram and the synthetic cluster sequences computed from stellar 
evolutionary models. In the conventional isochrone fitting procedure, the age is esti- 
mated in a subjective way. Here, we calculate the degree of likelihood by applying 
a modern statistical estimator, the Saha estimator W, and the interval of confidence 
from x 2 statistics. We apply this approach to sets of three different evolutionary models. 
Each of these sets consists of different chemical abundances, ages, input physics, etc. 
Based on our approach, we estimate the age of NGC 6397, M92 and M3. With a con- 
fidence level of 99%, we find that the best estimate of the age is f4.0 Gyr within the 
range of f 3.8 to 14.4 Gyr for NGC 6397, 14.75 Gyr within the range of 14.50 to 15.40 
Gyr for M92, and 16.0 Gyr within the range of 15.9 to 16.3 Gyr for M3. 



1 Introduction 

Since the Globular Clusters (GCs) contain the oldest objects for which age esti- 
mates are available, they have been recognized as being of key importance for de- 
riving the age of the Universe. In some conventional procedures for determining 
the age of GCs, like the isochrone fitting, the stellar evolutionary model has 
been selected in a subjective way. There, the observed position of the stars in 
the Colour-Magnitude Diagram (CMD) is compared with the position calculated 
from theoretical models. The model which minimizes the discrepancies between 
the observed and calculated sequences gives the estimated age of the cluster. 
Motivated by these ideas, and due to the refinement of theoretical stellar models 
and to a growing number of observations, we present a statistical approach JlT[ 
that allows to estimate the age in a more objective way by using the Saha statis- 
tics W Q and x 2 0- As an example, we subsequently derive the age of NGC 
6397, M92 and M3. 



2 Observational data and stellar evolutionary models 

The observational data consist of V and I photometry of three GCs of the Galaxy: 
NGC 6397 D [§, M92 and M3 [pj. We removed stars with excessive photometric 
errors (3er in the colour) and we selected only stars out to a radius of r<140' for 
M92 and 170' for M3. After the procedure of rejecting stars, the sample of M92 
comprises 4846 stars and the sample of M3 10333 stars (Fig. II]). 
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Fig. 1. CMDs of NGC 6397, M92 and M3, after the procedure of removing stars and 
considering the adopted parameters given in jip[| 

We selected the models computed by Q and (t]] (Padova isochrones), tracks 
developed by || (Yale isochrones) and models by || Q (Pisa isochrones). We 
chose them as they were published in the theoretical plane (HRD). Every model 
is computed for different input data (metallicities, ages, opacities, stellar atmos- 
pheres, etc.), resulting in differences between the isochrones. On the other hand, 
every model uses different forms of conversion between the theoretical and the 
observational plane. Because this lack of uniformity may influence the estimation 
of the age, we limited this additional uncertainty by using the same transforma- 
tion to each set of isochrones: tabulations of bolometric corrections and colours 
determinated from the SED extracted from the spectral library of Kurucz [0] . 



3 The method and the age of NGC 6397, M92 and M3 



Let S and M be the number of observed stars and synthetic stars generated from 
a particular stellar model in CMDs of a GC. If both CMDs are binned in a grid 
of B bins, the probability of having gathered the data under a particular model 
is the likelihood function, which is given in |l|. By mi we denote the number of 
synthetic stars placed on the i th bin and by Si the number of observed stars 
inside of the i th bin. 

Since there is not a unique configuration of a synthetic CMD of a GC, we 
generate 1000 different synthetic configurations of the CMD originating from 
every isochrone of age ij. Subsequently we draw the observational CMD and the 
synthetic CMDs within a grid of B cells. We then count the number of stars 
in each cell for every CMD. In the next step we first calculate log W between 
two synthetic CMDs (generated from the same isochrone of age ti). Secondly, 
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Fig. 2. Example of the distribution of values of log Wfor different ages of the Padova 
isochrones for a sample of 1187 stars in NGC 6397 (see text for details) 

we estimate log W for a new couple of different synthetic samples and continue 
determining the value of log W between different samples until 500 couples are 
completed (model-model comparison). The values computed of log W give a 
distribution of log W, which is represented as the non-shaded region in the 
histogram (Fig. ||). 

We calculate the model-model distribution to compare it with the distri- 
bution of frequencies of log W determinated from comparing observational and 
synthetic samples (data-model comparison). The data-model distributions are 
constructed by computing log Wiot couples between 500 synthetic samples (from 
the same isochrone of age ti) and the observational data. In this way 500 values 
of log W are computed. The distribution of log W is shown as the shaded area 
in the histogram of Fig. Q. 

We estimate the interval of confidence by using the \ 2 statistical test, as 
suggested by |9j, between the observational and synthetic data and for every 
evolutionary model. We chose the lowest value of % 2 , Xmin- The quantity in- 
creases when we move away from the best fit. The values of the confidence levels 
for a parameter and its respective probabilities (in percent) for data distributed 
in form of a Gaussian curve were taken from the statistical tables of The 
estimate of the age is given by the age of the isochrone that corresponds to 
the minimal distance obtained between the Gaussian median fits to the model- 
model and data-model distributions (14.75 Gyr in Fig. ||). The full details of the 
approach are presented in (hJ. 

Sl(M + B - 1)! 

prob(data | model) = {M + S + B _U [ ' w W 
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where W is given by 

i=\ 11 i=l v 1 ' 



Table 1. Summary of the results of the best estimate of the age of the samples con- 
sidered in this work, for every evolutionary model 



Evolutionary 


Z of the 






Best estimate Interval 


of 99% 


Model 


Isochronc 






of t [Gyr] 


of confidence 


NGC 6397 


1187 stars 


B=14400 










Padua 


0.0004 


1352.81 


2.82 


14.50 


[14.3 - 


15.1] 


Yale 


0.0004 


1089.04 


2.34 


14.00 


[13.8 - 


14.3] 


NGC 6397 


373 stars 


B=8400 










Padua 


0.0004 


232.84 


1.41 


14.25 


[13.7 - 


15.6] 


Yale 


0.0004 


206.80 


1.25 


14.00 


[13.3 - 


15.1] 


Pisa 


0.0002 


249,95 


1.41 


13.00 


[11.9 - 


14.2] 


M92 


4846 stars 


B=12000 










Padua 


0.0001 


2221.63 


1.31 


14.75 


[14.5 - 


15.4] 


Yale 


0.0002 


2455.26 


1.42 


15.00 


[14.6 - 


15.7] 


M92 


1482 stars 


B=12000 










Padua 


0.0001 


678.79 


1.73 


15.00 


[14.8 - 


15.5] 


Yale 


0.0002 


769.67 


1.91 


15.00 


[14.7 - 


15.5] 


Pisa 


0.0002 


832.86 


2.06 


12.00 


[11.8 - 


12.1] 


M3 


10333 stars 


B=14400 










Padua 


0.0004 


3558.38 


1.59 


15.75 


[15.7 - 


15.9] 


Yale 


0.0004 


3037.01 


1.37 


16.00 


[15.9 - 


16.3] 


M3 


4929 stars 


B=14400 










Padua 


0.0004 


1038.53 


2.14 


16.25 


[16.2 - 


16.6] 


Yale 


0.0004 


757.20 


1.68 


16.00 


[15.9 - 


16.3] 


Pisa 


0.0002 


1120.65 


2.45 


14.00 


[13.9 - 


14.2] 
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